Search for: All records

We present a new method to estimate the boundary of extended sources in high-energy photon lists and to quantify the uncertainty in the boundary. This method extends the graphed seeded region growing method developed by M. Fan et al. Here, we describe how an unambiguous boundary of a centrally concentrated astronomical source may be defined by first spatially segmenting the photon list, then forcibly merging the segments until only two segments—an extended source and its background—remain, and finally constructing a boundary as the connected outer edges of the Voronoi tessellation of the photons included in the source segment. The resulting boundary is then modeled using Fourier descriptors to generate a smooth curve, and this curve is bootstrapped to generate uncertainties. We apply the method to photon event lists obtained during the observations of galaxies NGC 2300 and Arp 299. We demonstrate how the derived extent and enclosed flux of NGC 2300 obtained with Chandra and XMM-Newton are comparable. We also show how complex internal structure, as in the case of Arp 299, may be subsumed to construct a compact boundary of the object. 

The production of complex astronomical data is accelerating, especially with newer telescopes producing ever more large-scale surveys. The increased quantity, complexity, and variety of astronomical data demand a parallel increase in skill and sophistication in developing, deciding, and deploying statistical methods. Understanding limitations and appreciating nuances in statistical and machine learning methods and the reasoning behind them is essential for improving data-analytic proficiency and acumen. Aiming to facilitate such improvement in astronomy, we delineate cautionary tales in statistics via six maxims, with examples drawn from the astronomical literature. Inspired by the significant quality improvement in business and manufacturing processes by the routine adoption of Six Sigma, we hope the routine reflection on these Six Maxims will improve the quality of both data analysis and scientific findings in astronomy. 

Data from high-energy observations are usually obtained as lists of photon events. A common analysis task for such data is to identify whether diffuse emission exists, and to estimate its surface brightness, even in the presence of point sources that may be superposed. We have developed a novel nonparametric event list segmentation algorithm to divide up the field of view into distinct emission components. We use photon location data directly, without binning them into an image. We first construct a graph from the Voronoi tessellation of the observed photon locations and then grow segments using a new adaptation of seeded region growing that we call Seeded Region Growing on Graph, after which the overall method is named SRGonG. Starting with a set of seed locations, this results in an oversegmented data set, which SRGonG then coalesces using a greedy algorithm where adjacent segments are merged to minimize a model comparison statistic; we use the Bayesian Information Criterion. Using SRGonG we are able to identify point-like and diffuse extended sources in the data with equal facility. We validate SRGonG using simulations, demonstrating that it is capable of discerning irregularly shaped low-surface-brightness emission structures as well as point-like sources with strengths comparable to that seen in typical X-ray data. We demonstrate SRGonG's use on the Chandra data of the Antennae galaxies and show that it segments the complex structures appropriately.